BCS-BEC crossover in a relativistic superfluid and its significance to quark matter 
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The character change of a superfluid state due to the variation of the attractive force is inves- 
tigated in the relativistic framework with a massive fermion. Two crossovers are found. One is 
a crossover from the usual BCS state to the Bose-Einstein condensation (BEC) of bound fermion 
pairs. The other is from the BEC to the relativistic Bose-Einstein condensation (RBEC) of nearly 
massless bound pairs where antiparticles as well as particles dominate the thermodynamics. Possible 
realization of the BEC and RBEC states in the quark matter is also pointed out. 
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Recently, new superfluid states in the ultracold gas of 
fcrmionic alkali atoms ( 40 K, 6 Li) were realized Using 
the Feshbach resonance, the long-standing idea of the 
crossover from the BCS state to the Bose-Einstein con- 
densation (BEC) [3- El- [Hi has been extensively examined. 
The basic concept of the BCS-BEC crossover is as follows: 
As long as the attractive interaction between fermions is 
weak, the system exhibits the superfluidity characterized 
by the energy gap in the BCS mechanism. On the other 
hand, if the attractive interaction is strong enough, the 
fermions first form bound molecules (bosons) . Then they 
start to condense into the bosonic zero-mode at some 
critical temperature. These two situations are smoothly 
connected without the phase transition. 

The possible realization of the BCS-BEC crossover 
in various systems has been theoretically investigated. 
These include the liquid 3 Hej3| , the trapped alkali atoms 
||| , and the nuclear matter U One of the most striking 
features of the crossover is that the critical temperature 
in the BEC region is independent of the coupling for the 
attraction between fermions. This is because the increase 
of the coupling only affects the internal structure of the 
bosons, while the critical temperature is determined by 
the boson's kinetic energy. Thus, the critical tempera- 
ture reaches a ceiling for the large coupling as long as 
the binding effect on the boson mass can be neglected. 
Even in the nuclear matter where the interaction is rela- 
tively strong, the binding energy of the deuteron is much 
smaller than the nucleon mass. This fact allows us to 
work within a nonrelativistic framework for describing 
such a crossover. 

It is interesting to ask how the situation changes in 
relativistic systems where the binding effect can not be 
neglected. The color superconducting phase in the dense 
quark matter 0,13 and the pion superfluid phase at finite 
isospin density |9| would be examples. In this article, 
we will show that there could be two crossovers in the 
relativistic superfluids. One is the ordinary BCS-BEC 



crossover, where the critical temperature in the BEC re- 
gion would not plateau because of the relativistic effect. 
The other is from the BEC state to the novel state, the 
relativistic BEC (RBEC), where the critical temperature 
increases to the order of the Fermi energy. 

In order to explore the BCS-BEC and BEC-RBEC 
crossovers in the relativistic system, we start with the 
following contact four-Fcrmi interaction model: 



£[ip, tp]=ip (i@ - m + 70^) ip 

+ G (i^ 15 CiJj*) ■ (# T 75^) 
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Here if> = ifi(t, x) is a Dirac fermion field having a spinor 
index implicitly. The fermion mass and chemical poten- 
tial are denoted by m and ji. C = 17072 is the charge 
conjugation matrix and G is a coupling constant for the 
attraction in the J p = + channel. Qualitative results 
shown below are not modified even when fermions have 
internal degrees of freedom other than spin. Therefore, 
we will make our analysis without them for simplicity. 
The partition function can be written as 

Z = J VipVip exp (~ ' dr J dxC^,i)]j , (2) 

where £e is the Lagrangian density in the Euclidean 
space. Introducing Hubbard-Stratonovich fields A(r, a;) 
for i%p T ^ 5 Cip and A*(r, x) for iip^^Ctp* and integrating 
out the fermion fields lead to 



Z = Z Q J VAVA* exp (-SW [A, A*]) . 



(3) 



Here Zq = e _/3sl °( AI ' T ) is the free fermion part of the par- 
tition function, while S ff[A,A*] is the effective action 
for the bosonic fields. In order to include the effect of 
the fluctuation, we evaluate the functional integral in the 
Gaussian approximation, whose validity will be discussed 
later. Expansion of the effective action up to the second 
order in A results in 
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where A(iui n ,p) is the momentum representation of the 
pair field with w„ = lirnT being the Matsubara fre- 
quency. x(ik>n,p) is a pair susceptibility at the one loop 
level [l(J]. The critical temperature for the superfluidity, 
T c , is given by the solution of the equation: 



^-x(o,o) 



= 0. 



T=T C 



(5) 



This is nothing but the Thouless criterion which states 
that the pair fluctuation becomes tachyonic at low mo- 
mentum because of 1/G — x(0, 0)\t<t c < 0. This is the 
signal of the BCS instability to the formation of Cooper 
pairs. 

The integration over A, A* in Eq. leads to the ther- 
modynamic potential in the Gaussian approximation: 
T) = f2 (A<, T) + fifl uc (/i, T) with 



^fluc(/*t, t) = t 



dp 



log 



1 

G 



X{iu n ,p) 



(6) 



Following Nozieres and Schmitt-Rink Q, f2g uc can be 
written in terms of a phase shift 6(u>,p) defined by 
S(u!,p) = — arg [1/G — x( w + *0,p)]. By differentiat- 
ing the thermodynamic potential with /i, we obtain the 
fcrmion number density as follows: 



iVtotal 



dp 



{f F (E p -»)-f F (E p + (i)} 



(2tt) 3 

f dp [°° dcu 85 



(7) 



with /f(w) = l/(exp[w/T] + 1) being the Fermi distri- 
bution function and /b(^) = sign(w)/ (exp[|cj| /T] — 1) 
being the Bose distribution function. 1 The first term 
which we denote by -/Vmf = Np — Np represents the con- 
tribution of fermions and antifermions at the mean field 
level and the second one which we denote by -/Vfl uc rep- 
resents the contribution of pair fluctuations. Instead of 
•Ntotai, we w iH sometimes use the Fermi momentum pp, 
which is defined by AT tota j = pp/3n 2 . 

If the attraction is strong enough, bound states ap- 
pear and we can extract the bound boson (antiboson) 
contribution Nq (A/g) from iVfl uc 0,O|. By picking up 
the bound state poles in d6(uj,p)/dn in the w-integral of 
Eq. J7J, we obtain 
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/b (wb(p) - 2/x) (8) 



h{u^{p) + 2 i jL). (9) 



Here, u>b(p) and — wg(p) are the solutions of 1/G — — 
2/x,p) = and correspond to the energy of the boson 
and the antiboson, respectively. Then the remaining part 
^unstable = Nfiuc _ ( Nb _ N ^ can be interpreted as the 

contribution of unstable off-shell bosons. 

In numerical calculations, a momentum cutoff A is 
introduced in order to regularize the ultraviolet diver- 
gence and all the dimensionful quantities are scaled by 
A. We take a characteristic parameter set (m/A = 0.2 
and pf/A = 0.1) so that we can analyze the effect of 
relativity. We confirmed that the variation of p F does 
not change our qualitative arguments below. Also, how 
the variation of m affects our results will be discussed 
later. Figures Ufa) and ^b) show numerical results of 
the critical temperature T c and the chemical potential fi 
as functions of G with the total number density A^totai 
fixed, which are obtained by solving Eqs. JSJ and Q si- 
multaneously. The ratios of the fermion and stable boson 
densities to ATtotai at T = T c are also plotted in Fig.^c). 
Based on these three figures, we will argue that there are 
three physically distinct regions; the weak, intermediate, 
and strong coupling regions. The superfluid states real- 
ized in the three regions will be interpreted as the BCS, 
BEC, and relativistic BEC phases, respectively. 

In the weak coupling region G/Go < 0.86, T c in- 
creases exponentially as is well-known in the weak cou- 
pling BCS theory. Its behavior is well described by 
the mean field approximation (the left thin solid line in 
Fig. Ola)), /i in this region is almost equal to the Fermi 
energy E-p — \J m 2 + pp . Accordingly, the fermion den- 
sity Np dominates the total density. From these facts, the 
superfluid state realized in this region can be regarded as 
the BCS state. 

In the intermediate coupling region 0.86 < G/Go < 
1.07, T c increases much slowly and fi decreases 
monotonously. Once /i becomes smaller than m, sta- 
ble bosons with the mass Mb(T c ) = 2/i appear and they 
dominate the total density. 2 The critical temperature for 
the ideal Bose gas is approximately given by 



rriNR 

J BEC 



2tt 
Mb 



1 2/3 



2C(3/2) 



(10) 



This nonrelativistic formula for the BEC critical temper- 
ature with the boson mass Mb(T c ) = 2/i is examined by 
the dotted line in Fig.^a), which well approximates T c 
in the intermediate coupling region. Therefore, we can 
interpret the superfluid state realized in this region is in 
the BEC phase. The slowly increasing T c in the BEC 
phase is in contrast to the result of the nonrelativistic 
calculation where the critical temperature approaches to 
a constant value |ll|. In the nonrelativistic framework, 
the change of the boson mass due to the binding can be 



2 The apparent singularity in the stable boson density in Fig.HTc') 
1 We use this Bose distribution function so that the thermody- does not mean a phase transition. The total boson density Nr uc 

namic potential is symmetric under fi — » — /i. is a smooth and positive function of the coupling. 
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neglected by definition. Thus, the critical temperature in 
the BEC phase is given by with Mb = 2m, which 

is indicated by the arrow in Fig.^a), and is independent 
of the coupling. In our rclativistic framework, however, 
the boson mass Mb(T c ) = 2/x can become smaller, and 
consequently, T c becomes larger as one increases the cou- 
pling. 

We remark that the criterion for the BEC state, \i < 
to, is model independent. At T = T c of the BEC, the 
chemical potential for the boson /ib = 2/x should be equal 
to its mass Mb- On the other hand, Mb must be less 
than 2to for the binding. Therefore, we have \i < to in 
the BEC region. 

Let us discuss the strong coupling region G/Gq > 1.07, 
where T c rapidly increases and the nonrelativistic formula 
for the BEC critical temperature breaks down. Because n 
is smaller than T c in this region, antiparticles can be eas- 
ily excited. As is shown in Fig.^c), the antifermion and 
antiboson densities (iVp,iVg) grow rapidly. At the same 
time, the fermion and stable boson densities (Np,Nb) 
increase so that the total number density is unchanged. 
T c in this region can be approximated by the ideal BEC 
critical temperature in the relativistic limit |l2^: 



rriRL _ 

j bec — 



'3(AT B -iV B ) 



2Mb 



(11) 



We note this approximate formula slightly deviates from 
T c particularly in the large coupling region. This is be- 
cause a large number of fermions are accompanied there, 
which is favorable in terms of entropy. We refer to the 
boson condensed phase with antiparticles in the strong 
coupling region as the rclativistic BEC (RBEC) phase. 

We can also understand the rapid increase of T c in 
terms of the decreasing mean interparticle distance d p . 
Let us estimate the critical temperature in the RBEC 
phase by comparing the thermal de Broglie wavelength 
7r/(-\/3T) to d p . We estimate d p by Np since fermions 
give a dominant contribution to the density in the strong 
coupling region. Thus, we have T ~ ttN^/Vs, which 
agrees with T c well (see the right thin solid line in 
Fig. da)). 

An essential difference between the BCS and (R)BEC 
phases is that the stable bosons are present above T c in 
the (R)BEC region. As the temperature is increased, 
the binding energy of the stable boson decreases. The 
bound state pole eventually disappears at a certain tem- 
perature, which we call a dissociation temperature Tdi ss - 
Tdi S s as a function of the coupling is shown in Fig. ^a) 
by the dashed line. Tdi ss line appears from the point 
G/Gq ~ 0.92 where the fermion pairs start to form 
the bound bosons. They get bound deeper with the in- 
creasing coupling, and as a consequence, Tdi ss increases 
monotonically. Tdi ss line separates the normal phase into 
two regions; a normal phase without stable bosons for 
T > Tdiss, and a preformed boson phase with stable 
bosons for T c < T < Tdiss- The preformed boson phase 
in the intermediate and strong coupling regions may pro- 
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FIG. 1: (a) Critical temperature T c normalized by the Fermi 
energy Ef — U m? + p F (thick solid line) as a function of 
the coupling G/Go- Go is defined by Go = 7r 2 /A 2 . For other 
lines, see the text, (b) Chemical potential fi/Ep as a function 
of the coupling G/Gq. The dashed line represents the level 
where fi = m; m/Ev ~ 0.89 in the present case. The dotted 
line corresponds to one half of the bound boson's mass in the 
vacuum, (c) The ratios of Nf,N-p,Nb and iVg to the fixed 
A/total as functions of G/Gq. N M f /Wtotai and AW/j'Vtotai are 
also plotted. The line for TVmf is behind that for Nf for 
G/Go < 1.07 because of the absence of antifermions. 



vide us with new insight into the precursory phenomenon 
above T c [13 . 

We have discussed the character change of the super- 
fluid state with a specific set of to and p-p. Let us now dis- 
cuss the m-dependence of the crossover boundaries with 
keeping p-p fixed. The crossover boundary from the BCS 
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to the BEC is characterized by the point where the bound 
states are formed in the medium. On the other hand, the 
crossover boundary from the BEC to the RBEC is char- 
acterized by the point where the boson mass 2fj, becomes 
smaller than T c . These two points are well approximated 
by the coupling G\ where the boson is formed and G2 
where it becomes massless in the vacuum with ji = T = 
(see the dotted line in Fig^b) which is indistinguishable 
from the solid line). We can show that Gi(2) increases 
(decreases) as decreasing the fermion mass m. It means 
that it becomes hard to bind two fermions due to the 
larger kinetic energy for a smaller m, while less attrac- 
tion will be needed to cancel 2m by the binding. Accord- 
ingly, the BEC region shrinks with decreasing m, while 
the RBEC dominates the larger region in the coupling 
space. In the ultrarelativistic limit m — > 0, the BEC re- 
gion will disappear because fermions could no longer be 
bound. Even in this case, we still have a superfluid phase 
with 2/i < T c for the large coupling which is smoothly 
connected with the RBEC phase at m > 0. 

In summary, we have discussed two crossovers in the 
relativistic four-Fermi model with the massive fermion: 
One is the crossover from the usual BCS to the BEC 
of bound fermion pairs and the other is that from the 
BEC to the RBEC of nearly massless bound pairs. In 
order to avoid the cutoff artifacts, we have checked each 
of Np, N-q, Np and Ag does not exceed (A/27r) 3 within 
the coupling range shown here. 

Since we have employed the Gaussian approximation 
which corresponds to the resummation of ring diagrams 
into the thermodynamic potential, all the interactions 
among bosons are neglected in our analysis. However, 
it is shown in the nonrelativistic framework that the 2- 
body interaction between bosons becomes smaller with 
increasing the attraction between fermions The 
multibody scatterings among bosons are also negligible 
in a dilute gas. Thus, our approximation is valid except 
for the vicinity of the BCS-BEC crossover boundary and 
the very dense RBEC limit. Going beyond the Gaussian 
approximation so as to take into account the interactions 
among bosons in our relativistic framework is an impor- 
tant future issue. Other approaches to the BCS-BEC 
crossover also may be useful [l3| . 

Finally, we make some speculative remarks on the rel- 
evance of the phases discussed above to QCD. The BCS- 
BEC crossover, which takes place for > T, may be 
realized in the cold dense quark matter || . The fermion 
mass m in Eq. in this case should be interpreted as 
the current or dynamical quark masses. Also, it is an in- 
teresting future problem to generalize our model by tak- 
ing into account the plasmino mass m p i ~ Sv^M 2 + tt 2 T 2 
with g being the QCD coupling constant. The plasmino 



mass can play a role of the chiral invariant mass consti- 
tuting the boson mass. Further study with the plasmino 
effective action will give us more insight into the re- 
alistic BCS-BEC crossover in the quark matter. In fact, 
the BEC criterion /i < m p \ leads to g > 1 for /1 ^> T, 
which corresponds to the density relevant to the center 
of compact stars. 

The BEC-RBEC crossover, which takes place for T 
fi ~ 0, will be relevant to the quark-gluon plasma just 
above the deconfinement transition. Possibility of having 
not only qq bound states [15j but also qq bound states 
in the deconfined phase has close relevance to the 
bound bosons in the RBEC state in this article. The 
plasmino mass m p i ~ gT again will play a crucial role to 
have bound bosons in the realistic situation. 

The above discussions suggest that there is a band 
of superfluid phases (BCS-BEC-RBEC) between the 
hadronic phase and the quark-gluon plasma phase in the 
QCD phase diagram. Also, the preformed boson phase 
may exist between the (R)BEC phase and the quark- 
gluon plasma phase. The transport properties in these 
phases are also of great interest. In the nonrelativistic 
framework, it is shown that the pair fluctuation is dis- 
sipative in the BCS region, while it propagates without 
viscous damping in the BEC region The detailed 

analysis including all the relevant hydrodynamic modes 
may provide a picture for almost the perfect fluid aspect 
of quark-gluon plasma. 

The generalization of our work so as to allow for the 
qq condensation will be essential to see the realization of 
the BEC and RBEC states in QCD. It is known that the 
one-gluon exchange generates the attraction in the scalar 
qq channel whose strength is 2 times larger than that in 
the qq channel. Even though the qq condensation still 
has a kinematic advantage due to the existence of the 
Fermi surface, the dominant qq attraction may wash out 
the (R)BEC phase leading to the large qq condensation. 
Whether the (R)BEC phase survives in the QCD phase 
diagram should be settled after taking into account the 
possibility of the qq condensation in our analysis. 
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